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Abstract 

In this paper we studied the double scaling limit of a random unitary matrix 
ensemble near a singular point where a new cut is emerging from the support of the 
equilibrium measure. We obtained the asymptotic of the correlation kernel by using 
the Riemann-Hilbert approach. We have shown that the kernel near the critical point 
is given by the correlation kernel of a random unitary matrix ensemble with weight 
e -x 11 _ ^his provides a rigorous proof of the previous results in 



1 Introduction 

In this paper we studied a double scaling limit of the unitary random matrix model with 
the probability distribution 

Z^ N exp(-Ntr(V(M)))dM, Z n , N = j exp (-Mr( V(M)))dM (1.1) 

J Tin 

defined on the space 7i n of Hermitian n x n matrices M, where V is real analytic and 
satisfies 

V(x) 

lim i — / 2 , 7x = +oo- 
x^±oo log(ar + 1) 

The eigenvalues x 1; . . . , x n of the matrices in this ensemble is distributed according to the 
probability distribution (See, e.g. [25], [T2] ) 

V {n ' N \x u . . .,x n )d n x = Z-^ N ^Un-i) Y[( Xj - Xk f dxi . ..dx n , (1.2) 

j<k 

where Z n ^ is the normalization constant. 

A particular important object is the m-point correlation function lZm' N \xi, ■ ■ ■ ,x m ) 

TZ^ N \x 1 ,...,x m )= nl f ■■■ [ P^ N \x 1 ,...,x n )dx m+1 ...dx n . (1.3) 
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The correlation function is a very useful quantity in the calculation of probabilities. In 
fact the 1-point correlation function TZ^ 1 (x) gives the probability density of finding an 
eigenvalue at the point x. (Note that, however, the m-point correlation function TZm is 
not a probability density in general.) 

A well-known result concerning the m-point correlation function is that it admits a 
determinantal expression with a kernel constructed from orthogonal polynomials. (See e.g. 

mm 

To be precise, let ir n (x) be the degree n monic orthogonal polynomials with weight 
e -NV(x) Qn R> [29] 

n n {x)Tx m (x)e~ NV{x) dx = h n 5 nm . (1.4) 
Let us construct the correlation kernel by 

n(v(x)+v(x')) sr^ 7r j( x ) 7r i( a: 



hi 

3=0 3 

By the Christoffel-Darboux formula, this kernel can be expressed in terms of the two 
orthogonal polynomials n n (x) and 7r n _i(x) instead of the whole sum: 

K nN (x X') = c -hN{V(x)+V(x')) 7l n(x)n n ^ 1 (x') - TT n (x')TT n - 1 (x) ^ 

h n -i(x-x') 

Then the m-point correlation function can be written as the determinant of the kernel ( 11 .51) 
P2], [25], [27] 

Km N) ( x ii ■ ■ ■ i x m) = det (K^ixjiXk))^^ 

In the limit n,N oo, ~ 1, the 1-point correlation function ll±' N \x) of the 
ensemble ( II. ip is asymptotic to the equilibrium measure p(x) [12], [22], [28J : 

lim TZ^' N \x) = p(x), 

where the p(x)dx = dp, min (x) is the density of the unique measure p m in(x) that minimizes 
the energy 

1(h) = — / / log \x — y\dji(x)dp,(y) + / V(x)dfi(x) 



among all Borel probability measures p on M. The fact that p m i n (x) admits a probability 
density follows from the assumption that V(x) is real and analytic [Hj. Moreover, it was 
shown in p3] that for real and analytic V(x), the equilibrium measure is supported on a 
finite union of intervals. 

The following conditions are satisfied by the equilibrium density p(x) [12], [21] 



2 / log \x — s\p(s)ds — V(x) = I, x G Supp(p(x)), 

'.. " (1.6) 

2 / log \x - s\p(s)ds - V(x) <l, x G R/Supp(p(x)). 
Jr 
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For a generic potential V(x), the inequality in ( 11 .61) is satisfied strictly 

2 / log \x — s\p(s)ds — V(x) > I, x G R/Supp(p(x)). 
Jr 

However, for some special potential V(x), this inequality may not be strict and we may 
have 

2 / log \x — s\p(s)ds — V(x) = I, x = x* 
Jr 

at some point x* ^ Supp(p(s)). In this case, if we change the potential slightly, a new 
interval may emerge from the support of the equilibrium measure. This is the 'birth of 
new cut' critical limit that we are going to consider in this paper. 

According to [TSj, the 'birth of new cut' critical limit is a singularity of type I for the 
random matrix model (11.11) . Other singular cases include: 

1. Type II singularity is the case where the equilibrium density vanishes at a point x* 
inside the support. 

2. Type III singularity is the case where the equilibrium density vanishes faster than 
a square-root at an edge point x* of the support. (Generically it vanishes like a 
square-root at the edge) 

The asymptotic behavior of a random matrix ensemble near singular points has been 
studied extensively [2], @|, [5], [7], [8], [9], pig, pi], p3], [20], [30]. In these studies, one 
considers a one or multi-parameter family of potential Vt j (x) in which the singular point 
is achieved at tj = tj. One then studies the asymptotic behavior of the random matrix 
model (II. ip when tj is close to tj. The 'double scaling limit' is the study of the these 
asymptotic behavior when the differences between tj and tj are coupled with n and N. A 
remarkable feature is that in the double scaling limit, a universality can be observed. Upon 
a suitable scaling of the variables x and x', the asymptotic behavior of the kernel (11.51) near 
the critical point x* depends only on the type of singularity rather than the potential V(x) 
itself. 

In many cases, the behavior of the kernel in a double scaling limit is described by 
integrable hierarchies such as the Painleve equations. In the case of the type II singularity, 
[3], [H], [10] and [30] has shown that the kernel can be described by the Hastings-McLeod 
solution of the Painleve II equation in the double scaling limit. While for the type III 
singularity, the kernel can be described by the Painleve I transcendent [UJ, [13]. In [18] , 
the double scaling limit of the 'birth of new cut' was studied and the kernel was described 
by the orthogonal polynomials with weight e~ x v on the real axis. However, the formulae 
derived in [18] have not been rigorously proven and it is the purpose of this paper to provide 
a rigorous proof of these results. 
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V(x) 



i 



1.1 Statement of results 

We should now introduce some notations and state the results in this paper. 

In this paper, we should consider a one parameter family of potential V t (x) 
parametrised by t = jr. We should consider the double scaling limit of t — > 1 and n, 
iV — > oo such that 

lim ( ?L _ A = £/ . > 0, n > iV (1.7a) 

lim n fc f - = U- < 0, n < N, k € 1 - — , oo ) (1.7b) 

exist. In particular, for t < 1, we considered the regime where t — 1 is of order n~ k for any 
k greater than or equal to 1 — while the scaling for t > 1 is fixed. 

Let us now state the assumptions that are used in this study. Since the main point of 
this study is the treatment of the critical point x*, we will assume the followings: 

1. The support of the equilibrium density p(x) consists of one interval only, that is, 
the first equation of (jl.6p holds precisely on a single interval (a, b). Without lost of 
generality, we will assume that a = — 2 and 6 = 2. 

2. The equilibrium measure does not vanish at any interior point of (—2, 2). 

3. The point x* is the only point outside Supp(p)(x) where the inequality in ( 11. 6p is not 
strict and we assume that x* > 2. 

4. As pointed out in [T6j , the function 2 Llog \x — s\p(x)ds — V(x) — I vanishes to an 
even order at x*. We will assume that this order of vanishing is 2v. 

Let the equilibrium measure of V t (x) be p l (x) such that 

2 /loglx -.!✓(.)* -VJ(*)=fc, ^SuppM*)), 
Jr 

2 / log \x — s\p {s)ds - V t {x) <l t , x G R/Supp(p*(x)), 
Jr 

and denote by c x * the following 



'U 



n 



+ jj{s)\og\x*-s\ds^ >0. (1.9) 



It is known that both tp l (x) and the support of p t (x) are increasing with t [21], [13], [28] . 
[31] . In particular, for t < 1, the equilibrium measure is supported on one interval while 
for t slightly greater than 1, the equilibrium measure is supported on 2 intervals. 
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Let St be the support of p t (x). Then in [6], it was shown that the equilibrium measure 
d/j,t(x) = p t (x)dx satisfies the Buyarov-Rakhmanov equation 



to = - u Sr dr, (1.10) 
t Jo 

where uo$ T is the equilibrium measure of the set S T . Namely, it is the unique probability 
measure supported on S T that minimizes the logarithmic potential 



= J J — log \s — t\dfl(s)dfl(t) 



among all the Borel probability measures [i supported on S T . 
If S T consists of one interval only, then ujs t [x) i s given by 

uj St = — 7= = ==dx, x e (a T , b T ). 
7ry [b T — x)[x — a T ) 

In particular, we have, at t — 1 

hm = dx = w(x)dx. (1.11) 

t-,1 t-1 TTv/4^2 ^ ^ > 

The fact that w(x)dx is the equilibrium measure on the interval [—2, 2] means that 

f 2 S 

/ w(s) log \x — s\ds = -, i 6 [-2,2], 

(1-12) 

/ w(s) \og(x — s)ds = \ogx + 0(1), x — > oo 

for some constant <j. 

Let us defined a function 0(x) that is closely related to w(x)dx. 

0( x ) = i+/" w {s) log{x* - s)ds. (1.13) 
2 7-2 

In this paper, we will use an anzatz in [18J to construct an approximated equilibrium 
density p t {x) for t > 1 and use it to modify the Riemann-Hilbert problem of the orthogonal 
polynomials (jl.4p . 

We shall denote the correlation kernel for the random matrix model 

Z- 1 ! ,exp(-tr(M 2iy ))rfM, Z m>1/ = I exp (-tr M 2u )dM (1.14) 
by K^(x, x'). That is, 
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where vr^(x) is the degree m monic orthogonal polynomial on R with respect to the weight 
e~ x " and h v m is the corresponding normalization constant as in (11.4ft . 
We can now state our main result. 

Theorem 1.1. Let V{x) be real and analytic on R such that lim^-too = +oo. Let 

p(x) be the density of the equilibrium measure ofV(x) supported on the interval [—2,2]. 
Then 

where x* > 2 and Q(x) is real analytic on R with Q{x*) > 0. 

Let n, N — > oo siic/j £/ia£ ( [i. 7a| ) and ( i. 76| ) /jo/d and /ei it, w 6e i/ie following 



u = 2vcj){x*)U + 

(1.16) 



it 



2u<l>(x*)U + + \ 



where <fi(x*) is defined in M.ltft) and [x] is the greatest integer that is smaller than or equal 
to x. 

Let K n7 N be the correlation kernel ( 11. 5j) . i/ien /or u ^ N + |, i/ie /zmzt o/ the kernel is 
given by 

lim 1 K n>N lx* + Z -^,x* + ?-^)=K£{z,z'), w>JV, (1.17a) 

n,iv->oo (p(x*)n 2 " \ (p(x*)n 2 » (p(x*)n< 



\2u 



z ^ z' 



lim e c **K nN x* H -, x* + 



Z 2 »+U') 2V 1/1 1 

n < N. 



;i.l7b) 



87r \x* — f3 t x* — a ti 
where K^-(z,z f ) is defined in ( ti. 15\) and c x * is defined in ( 11.5)) and v 5 ^*) zs gu>en fry 



g(x*)v/(x*) 2 -4\ 2 ' 



2// 



The result shows that for u ^ N + |, the correlation kernel near x* for t > 1 is given 
by the correlation kernel of a finite random matrix ensemble (11.141) with size [u + -A . This 



confirms the results in |18j . When u goes pass a half integer, the size of the finite random 
matrix ensemble jumps by 1 and a non-trivial transition takes place. This is due to the 
non-uniform converges of ( I1.17al) in u when u is close to a half integer. When u is close 
to a half integer, error terms that depends on K£ ±l which are not seen in (11.1 7a[) become 
significant and start taking over the terms, which results in a jump when u goes pass 
a half integer. 



Birth of new cut double scaling limit 



7 



Note that (11.17bl) implies that the leading order term of the kernel at x* is e Cx * . This 
leading term tends to zero when n, N —>■ oo unless t — 1. This is not surprising as for t < 1, 
there is no eigenvalue near the point x* and the correlation kernel should be vanishing near 
x* in the limit. 

Remark 1.1. Claeys |?J/ has simultaneously and independently used the Riemann-Hilbert 
method to study the birth of new cut double scaling limit. In Claeys [1], the case when v = 1 
was studied and the Hermite polynomials was used to construct the asymptotic kernel. 
Despite the similarity of our work to [?J/, a very different treatment to the equilibrium 
measure was used in ^jj. In ^jj, the equilibrium measure with total mass 1 — 2^^<p(x*) 
was used to construct the 'g -function' for the Deift-Zhou steepest decent method when t > 1. 
Whereas in this paper, we approximated the equilibrium measure by solving the Buyarov- 
Rakhmanov equation ( ti. 10\) up to a certain order in t — 1. We then use this approximated 
measure to construct the 'g -function' for the Deift-Zhou steepest decent method. Also worth 
remarking is that in ^jj, the behavior of the kernel when u is close to a half integer was 
studied. 

This paper is organized as follows. In section [2] we will use the ansatz obtained in [18] 
to construct an approximated equilibrium density for t > 1. We then show that conditions 
of the type (jl.fip are satisfied for this approximated density outside some neighborhoods of 
the edge points and the critical point. We then study the error terms in these conditions. 

In section [3] we will apply the Deift-Zhou steepest decent method to the Riemann- 
Hilbert problem of the orthogonal polynomials ( 11. 4ft . We will use the approximated den- 
sity to construct a '(/-function' and use it to modify the Riemann-Hilbert problem. We 
then approximate this modified Riemmann-Hilbert problem by a Riemann-Hilbert prob- 
lem that can be solved explicitly and construct parametrices to solve this approximated 
Riemann-Hilbert problem. These parametrices then give us the asymptotics of the orthog- 
onal polynomials ( 11.41) . These asymptotics will then be used to derive the asymptotics of 
the kernel ( 11.51) in section HI 

2 Equilibrium measure 

We will now study the behavior of the equilibrium measure p t {x) (11. 8ft when t is close to 
1. Let t be a real parameter and let us define 

V t (x) = \v(x), t > 0. 

v 

Then V\{x) = V(x). We shall consider the case when t < 1 and t > 1 separately. For 
t > 1, we will replace the eigenvalues on the newborn interval by a point charge. Let the 
support of the equilibrium measure S t be 



S t = [a t ,b t ], t < 1 

S t = [a t ,b t ] U [ct,d t ], t>l. 



(2.1) 



8 



M. Y. Mo 



Let us define the function h l by 

h\x) = I log(x - s)dfn(s) (2.2) 



where the principal branch of the logarithm is taken in the above, 

log(x — s) — log |x — s\ + i arg(x — s) 

< arg(x — s) < 7r, s G R, 3x > 0, 
—Ti < arg(x — s) < 0, sGl, Sx, 0. 

The boundary values of /i*(x) on the real axis are then 

h±(x) = / log \x — s\dfi t (s) ± ni I dn t (s) 

JR J a t 

In particular, the function h l is analytic on C/[at, oo) and it satisfies the following 

h l + (x) + h l _{x) - V t {x) + l t = 0, x G [a t , 6 t ] U [q, d t ] 

h%(x) + - + l t < 0, x G E/ ([a t , 6t] U [cj, d t ] U {x*}) 

/" bt (2 3) 

h%(x) - ht(x) = 2ni dfi t (s), i£l 

/i*(x) = logx + 0(x _1 ) x — > oo 

In [T3], it was shown that for a real analytic potential V(x) on R, the equilibrium measure 
dfit{s) can be expressed in terms of the negative part of an analytic function qt(x). 

Theorem 2.1. fTlfl Let V(x) be real analytic in a neighborhood V of the real axis and let 
qt{x) be the following function 

V\x)V 1 f V'(x) - V\y) 



q t (x) = — ^ - - / —dfH(y), x G V. (2.4) 

Then the equilibrium measure has a density p t (x) which can be written as 

i\ v 

where q^{x) is the negative part of qt{x), that is, 

q t (x) = qt(x) + g t ~(x), g t + (x) > 0, q^(x) < 0. 
Moreover, we have the following 

- 6V - (2 - 5) 
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2.1 Approximated equilibrium measure for t > 1 

For t > 1, a new cut in the support of the equilibrium measure is emerging at x = x*. We 
would like to find an approximation to the equilibrium measure and study its properties. 

The Buyarov-Rakhmanov equation fll.lOp for the equilibrium measure is a nonlinear 
ODE which is difficult to solve. In [18], an ansatz was used to solve this differential 
equation up to some leading order terms in t — 1. As this ODE becomes singular at t = 1, 
it is difficult to prove rigorously that the solution in [18] does indeed give the equilibrium 
measure for t slightly greater than 1. 

Instead of showing that the solution obtained in [TH] gives the correct equilibrium 
measure for t > 1, we would use the ansatz in [TH] to construct an approximated density 
p f (x), together with a function h f (x) analogue to the function defined in (12.2ft . We 

will then show that this approximated density satisfies conditions of the type (11.61) up to 
a certain order in t — 1. 

First note that the function qt(x) defined in (12.41) has the following form at t — 1. 

y/tfx) = ~Q(x)(x - O^-V^^I, (2.6) 

where Q(x) is analytic in a neighborhood V of the real axis. 
We can now define a function qt(x) analogous to qt(x). 

Definition 2.1. Let 8t = t — 1 > 0. Then the function q l (x) is defined by 



VW) = y ^ a !! {X — \Q(x)H t (x)\l(x-x*y-4yi f-j^)' +V(x)6t 



(2.7) 



where a t and (3 t are, 



St St 
at = - 2+ (2 + x*) 2 - 1 g(-2)' fjt = 2 ~ (x*-2)2-iQ(2)' (2 ' 8) 

u>/w/e i?t(:r) is a monic polynomial defined by 

H t(x) = (* - g - XT" (-^) ' , (2.9) 

The function r)(x) is defined by 

Q(x)(x - X') 2 "- 1 Q{x){x - x*f»- 1 2_ 

V[X> 2Q(2)(2-x*) 2 »~ 1 (x-2) 2Q(-2)(2 + x*) 2 »- 1 (x + 2) x 2 -4' 1 } 

and the constant y is defined by 

^ Q(.r<) x /(.r) 2 - 1(2//)! 



v I | (2 . li; 



and 0(x*) defined in U.13\) . 
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Remark 2.1. The function n(x) is analytic in the neighborhood V of the real axis. 

We will now show that the density defined by the function \/qt(x) satisfies the Buyarov- 
Rakhmanov equation outside a fixed neighborhood of x* . 

Proposition 2.1. Let B s 5 be the set 

Bl = {x\ \x-s\<5} 

and let r\ = —2, r2 = 2 and r^ = x* . Then for sufficiently small 5t, there exist compact 
subset /C C V and S > independent on t, such that the function g'(x) satisfies 



V^l + °(^i)' -^/(y^U[-2,2]) (2.12) 
uniformly in V / fljj=i U [ — 2, 2] j , where St — t — 1. 



Proof. We will expand (12.71) in terms of 5t and — r^st- Let us first consider the product 



log St ■ 



H t (x)y (x — x*) 2 — Ay 2 ^— j^j^ Let S > be fixed. Then for small enough St, the 
following Taylor series expansion is valid outside of B'l . 

\\ {x - x*) 2 - Ay 2 f-r^-] " = y . t .^ j)l n .y j (x - x*)- 2j+1 (-r^r) " ( 2 - 13 ) 
V V log 5*/ ^j!j!(l-2j) y 1 ; v iog5ty 

Now from the Taylor series expansion of ^(x — x*) 2 — % 2 ^— , 



St \A 2 ^(2j)! _«_ af _i/ « 



logtft/ y ^ V. log«J* 

we see that (c.f. [IB] ) 

# t (x) = Pol ^(x - x*) 2 "- 1 ^(x - x*) 2 - V ( 5t 

where Pol(X) denotes the polynomial part of X. 
Therefore we have 



log St 



H t (x)\l (x - x*) 2 - Ay 2 (--f^y " = (x ~ x*) 2 "- 1 

CO , r, \ 1 + i 

+ E^ + ^-^- 1 (-r^) 'L,, (2.14) 



v-l 



log 5i 

(2p)! (2(j + i/-p)!) 

Q (J + u- p)\(j +u-p)\(l- 2(j +u-p)) 
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Then, for a small enough St, we have, for \x — x*\ > S, 
This means that, for x £ K,/B$ , we have 



H t {x)d{x - x*y - 4y* (-^) " = (* - iO 2 """ 1 + (^) > * £ 5 f ( 2 -15) 

Now let us look at the terms of order St. Again, for small enough St, the following Taylor 
series expansions are valid outside B\ U Bj 2 . 

, , ^ (-1V'(2j)! {E(-2)St\ j 

^ - 2j)4J \x-2j 

where the function is defined by 



(2.16) 



v ' (x- x*) 2v - 1 Q(x)' 
The identity (12. 16ft implies that, for small enough St, we have, for x £ K/ (B$ 2 U Bj), 



y/(x - a t )(x - A) = v^4 + ft ( - + ^MJ 1 ) + Omf)- (2-17) 

V 2\/x — 2 2v^ + 2 / 

Combining this with (12.151) and (12.101) . we see that, outside of Bf, the limit (12.121) is given 
by 



V^M-V^M ( v^T2S(2) v^+2S(2) 



t-1 \4E(x)Vx^2 4E(x)y/x=2 

(2.18) 

v^2H(-2) y^2H(-2) 1 \ ( St 



4E(x)Vx + 2 AE(x)Vx + 2 y/x 2 - 4 / V lo g^ 

which is just 



^ + °(l|«)< ^/Aj B? u[-2,2]) 



t-1 y/x* 

This gives the assertion of the proposition. □ 
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Let us now define the approximated equilibrium density to be 



p\x) 



— (V-<f (^)) > re e [a t , A] 



2tTX 



x £ [x* — a t , x* + a t ], af = 2y 



(2.19) 



St \ ' lv 



log St 



p\x) = 0, x G R/[a t , A] U [x* - a t ,x* + a t ], 



and let h l (x) be the following 



p (s) \og(x — s)ds + 



at 



X*+<Tt 



p (s) log(x — s)ds. 



(2.20) 



x*—at 



Then we have the following analogue of (11.111) for 



Proposition 2.2. For sufficiently small St, there exists S > such that the following is 
satisfied for h (x) 



~ f . . h(x) St 
hHx) = — ^ + — 
y ' t t 



wis) log(x — s)ds j + O 



St log(z + 2) 
log St 



xGC/U5?USupp(p*(x)), 



(2.21) 



where h(x) is the following 



h(x) 



p(s) \og(x — s)ds 



and w(s) is the equilibrium measure of the interval [—2,2] ( fi.iij) . 
Proof. Let us first divide the real axis in to different parts 

6 



U 

3=1 



where the 1L- are the following intervals, that is, 



Mi = [a t , -2 - 2 (2 + a t )] 



S S 
-2 + -,2 - - 
2' 2 



-2 - 2 (2 - at) , -2 + - 



2--,2-2(A-2) 



(2.22) 



^ 5 = [2 - 2 (A - 2) , A] , M 6 = [x - a t ,x + <r t ] 
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Let us now define T to be the line right above R3 

r = I x 



X = U + IE, U6 



-2 + S -,2- S - 
2' 2 



;2.23) 



Then we have 



p (s) \og(x — s)ds 



p(s) \og(x — s)ds 

3 

u>(s) log(x — s)ds 



tit 



log(x — s)ds, 



-q{s) 



TV 



r vrV4 



log(x — s)cis, 
= log(x — s)ds. 



(2.24) 



Let S > be such that the power series expansion of Q(x) and r^(x) around ±2 are valid 

inside Bf 2 . 

i 

First let us consider the integral on Mi. On Mi, the following power series expansions 
are valid. 

00 00 

q(*) = E ^o.- 2 )( s + 2 )"' = E ni* + 2 )' 

j=0 i=o 

00 

(\/& - s) = y/p t + 2^2x j (s + 2y, (2.25) 

00 1 / s + 2 V 
log(x - s) = log(x + 2) - E J > 

where the branch of log (a; + 2) is chosen to be the principal branch. 

It is not difficult to check that the coefficients in the above series remain finite as 5t — > 0. 
Moreover, from <\2. 15j) . we have 



HAx)\ (x - x 



*\2 



4y2 



St 



(x — X 



*\2u-l 



o 



St 



log St J \ — log St 

In particular, this means that on Mi the functions have the following estimates 

QO) = Q(-2) + 0(5t), r){s) = rj(-2) + O(St), 

%-s) = y//3 t + 2 + 0(6t), \og(x -s)= \og(x + 2) + 0(St). 

St 



, are Mi. 



H t (s)\ (s - x*) 2 - 4y 2 



St 
log St 



-2 + x 



*\2u-l 



o 



log St 
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Therefore the integral on Mi can be evaluated as 



p (s) \og(x — s)ds 



2S(-2)tt 
V / A + 21og(a; + 2)(S(-2))i 



y/s-a t ds(l + 0(5t)) 



(2.26) 



3tt 



(3St)* (1 + 0(5*)). 



Similarly, the following integrals for p(x) and the equilibrium measure on [—2, 2] are given 
by 



- 2 - 2<2+ "" ££) log( , _ s)ds = 2llo g (, + 2) 5 (-2)i , + 

2 I OI7T 

-2-2(2+a t ) 



W ( a )log(x- 8 ) da = log(x + 2) ^-^ + G 

t tlT 



(2.27) 



Therefore we have 

p*(s) log(a; — s) — 



- 2 - 2 < 2+ <*> / P ( S ) ft , 

U7(sj ) log (a; — sjas 



O ((5t)§log(x + 2)) . 



(2.28) 



Next let us consider the integral on K 2 - Since |x + 2| > 5, for s G M 2 , we can find constants 
independent on t and s such that 



\Q{s)\<M Q , \Js~+2\ < \Vs~=a~ t \ < M a , \rj(s)\<M v , 



Is — X 



* 1 2u- 



1 < M x *, \y/2=l\ <\y//3 t -s\< M p , 

S(2) 



2(2) 



s - 2 



<M 1; 



-.St 



\H t (s) - (s - x 



*\2u-l\ 



< M 



H 



< M 2 {St)\ 
St 



log 5t 



log | a; — s\ < M 3 log \x + 2|. 
Then, by using the the Taylor series expansion of y/s — at in fl 2 . X 6 1) . we see that 

1 



tp\s) - p(s) + St- 



ir 



V4^ 



log | a; — s\ 



< f £ V^+2(2j)! / |B(-2)|ft V 



Vs~+2(2j)\ f\E(-2)\6t 



-j\j\{2j - 1)4* V ^ + 2 



(2.29) 



„ v+ % /FT2(2j)! /|~(-2)|ftY / -5* \ \ , 
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for some positive constants E\, E2 and E%. One needs to be careful about the terms that 
contains negative power of x + 2 as they may become large in M2- If we integrate ( 12.29ft 
and consider the leading order term in St, we see that 



/ (tp t (s) — p(s) + Stw(s)) log \x — s\ds 
Jr? 



< E- 



St 



-log St 



log |x + 21. 



for some positive constant E. 
This implies that 



/ (tp\s) - p(s) + 5tw(s)) log \x - s\ds = O 

JRo 



We then see that 



St\og(x + 2) 
log 5t 

St \og(x + 2) 
log St 



(2.30) 



To compute the integral on K3, observe that for small enough St, the relation ( 12. 12ft holds 
uniformly on T. Therefore by (12.241) . the integral on R 3 is given by 



P 4 ( s ) - + —W(s) I lof> \.r - .s|J.s = () 



St log{x + 2) 
log St 



By applying the argument used for Ri and R2 to R4 and R5, we obtain 

'St \og(x + 2 



log St 



j =4,5. 



(2.3i; 



;2.32) 



Let us now consider the contribution from the interval [x — a t ,x + a t ]. From the power 
series expansions on R 6 , we have the following estimates, 



Q{s) = Q{x*) + 



St 



log St 



y^s - a t )(s - A) = V 7 ^*) 2 -4 + 
log(a; — sj = log(a; — x ) + (J 



r](s) = r](x*) + O 
St 



St 



log St 



1 



log5t 



log St 

1 - 



1 

2;/ 



;2.33) 



Therefore, the integral on satisfies the following estimate. 



r Q(^)y/(x*) 2 -41og(x-x^ 

/ P (s) log(x - sjrfs = 2 

Jr r ^ 



H t{s)\ a 2 t - (s - x*) 2 (is 1 + 



log 5t 



(2.34) 
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To evaluate the integral on the right, let us note that H t (x) can be written in the following 
form [T8] . 

St V~» „ ( , ^ ( St 



where P(s) is the following polynomial of degree 2v — 2, 

i:^/r^ i ;: - (2.36) 



Then by a change of variable 

C = (s - x* 



St N 2 * 



log <5t / 

in the integral on the right hand side of (12.341) . we have 

ft Q(x*)v/(x*) 2 -41og(x-x*) 

p (s) log(a; — s)as — 



log 5t 2tt 
' 2y . / ^ -\ . „ // St \- 

To evaluate this integral, we will use the following differential equation for P(£) in 



x 



/%K)(v^e) t 4> + o((^) 



(2.37) 



(2v - 2)P(0 - £P'(0 = - P(2y)) • 



Using this and integration by parts, we find that the integral in (I2.37P is given by 

-2y 



r m V^M=^^. (2-38) 

J-2V V 



Hence the integral (I2.37P is 



St y 2 P(2y)Q(x*)^(x*y - 41og(x - x* 

p (sj log (a; — sjds — 



log St v 
St 



x 1 + 



log St t 

Now by the use of induction, one can compute P{2y) easily [18J, 



u-l 



P(2y) = (2yr-^^4 



WK-k ..2,-2 (2f)! 



y 



k\k\ a 2{y-\)\iA 

3=0 ^ ' 
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This, together with the expression (12.111) for y implies that 

/ log(, - .)«(. = ~^") log(x - «•) (l + O ( (^) ) . (2.39) 



which is of order j^j. That is, 



Now by adding fl2~28l) . f[2~30D . pit . fl2~32l) . (E3Q]), we arrive at fl2~211 . □ 

Now from (11.12|) and (12.211) . we see that conditions of the type (12.31) are satisfied for 
h\x). 

Corollary 2.1. For sufficiently small 5t, there exist 5 > such that [x* — at, x* + <r t ] C Bf 
and that 

h* + (x) + hl(x) - ^ - ~ = (^) , x e [a t ,p t ]/ {Bf U Bl) 

» ^.(s) + ^» - ^ - < 0, x G R/ yj B% U [a*, A]j , (2-41) 

h\x) = (l + L t ^^JJlogx + 0(l), x^oo. 

where I is the constant I + (5t)q. The function v^(x) remains uniformly bounded in [ctt,Pt] 
as 5t — > 0, while the constant it remains finite in the limit. That is, if 

lim v+(x) = v (x), lim u = i. 

then v h (x) is uniformly bounded in [at, fit] an d i h is finite. 

Corollary 12.11 suggests that p'(x) is a good approximation to the actual equilibrium 
density p t (x). 

The following corollary follows immediately from the proof of proposition 12.21 and the 
fact that h(x*) - ^f 1 - | = 0. 

Corollary 2.2. Inside Bf , the following is satisfied. 

[* pHs) log(x - s)ds -^l = 5t F w(s) \og(x - s)ds + o( r^-] . (2.42) 

J at 1 J -2 VOgdtJ 

In particular, by using we see that the following is satisfied at x = x* . 

ft V(r*) I ( St \ 

f?(s) log(x* - s)ds - -±J- -- = 5t<P(x*) + O , (2.43) 

where <p(x*) is defined in U.13\) . 
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This corollary is essential for the construction of the local parametrix inside the neigh- 
borhood Bg* . (See section 13. 51) 



3 Riemann-Hilbert analysis 

A result by Fokas, Its and Kitaev [19] shows that the orthogonal polynomials (11.41) can 
be expressed in terms of a Riemann-Hilbert problem. In this section we will apply the 
Deift-Zhou steepest decent method to approximate this Riemann-Hilbert problem by a 
Riemann-Hilbert problem that is solvable explicitly. We will achieve this by using the 
approximated equilibrium measure constructed in section [21 We will modify the measure 
p t (x)dx by replacing the charges on [x* — at,x* + at] by a point charge. This then allows us 
to construct local parametrix near the critical point x* from orthogonal polynomials with 
weight e~ x " on the real axis. 



3.1 Riemann-Hilbert problem for the orthogonal polynomials 

One important property of the orthogonal polynomials (11.41) is that they can be represented 
solution to a Riemann-Hilbert problem [T9] . 
Consider the following Riemann-Hilbert problem for a matrix- valued function Y(x) = 
Y n)N (x). 

1. Y(x) is analytic on C/M. 

/i p -NV(x)\ 

2. Y + (x)=Y4x)[ Q 1 j, xeR (31) 

3. Y(x) = {l + 0(x- 1 ))^ x ° n y x^oo 

where Y + (x) and Y-(x) denotes the limiting values of Y(x) as it approaches the left and 
right hand sides of the real axis. This Riemann-Hilbert problem has the following unique 
solution. 



Y(x) 



n n{x) 2-wi Jr s-x 



where « n _i = —2nih n \ [12J. The correlation kernel (11.51) can be expressed in terms of the 
solution of the Riemann-Hilbert problem Y(x) via [H] 



KnA^v) = 2m{x _ y) (0 1) Y-\y)Y + (x) y (3.2) 

We shall apply the Deift-Zhou steepest decent method to the Riemann-Hilbert problem 
(13. ip to obtain the asymptotics for the orthogonal polynomials and the correlation kernel. 



V(x) + V(y) 
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3.2 Initial transformation of the Riemann-Hilbert problem 

We shall perform a series of transformation to the Riemann-Hilbert problem (13.1 1) and 
approximate it with a Riemann-Hilbert problem that can be solved explicitly. We will than 
use the solution of the final model Riemann-Hilbert problem to compute the asymptotics 
of the orthogonal polynomials and the correlation kernel. 

3.2.1 The ^-function 

To begin with, let us denote t by t = and rewrite the jump matrix in (13.1 ft as 

1 e -NV(x)\ /j e -nVt(x) 

1 ) = [o 1 

where Vt(x) = jV(x). We will now define a function g t (x) from the function h l (x) con- 
structed in section [2] to transform this Riemann-Hilbert problem. 
Let it* be the following 

u l = n I p t (s)ds 



where M 6 is defined in ( l2~22l) . 

For later convenience, let us denote by u 1 the non- negative integer closest to u l : 



7f 



u* > 

(3.3) 
= 0, u* < 0. 



From (12.391) . we see that if u* > 0, then 



i 



-* = -°i|s^( 1+0 ((^) Jj- (34) 

By inserting the scaling (II .7aj) into ( 13.41) . we see that u l is finite in this limit. 

As mentioned before, we would like to replace the charges in the interval [x* — a t , x* + a t ] 
by a point charge when t > 1. We should therefore define the ^-function to be the following. 



+ — log(z -x*), t > 1, ( 3 - 5 ) 
n 

g\x)= I p t (s) log(x — s)ds, t < 1, 

J at 
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where in the above equations, we have extended the definitions of the end points at and 
Pt (12.81) to include the the values of t that are less than or equal to 1: 

St St 
at= ~ + (2 + x*) 2 ^g(-2)' Pt ~ ~ (x*-2) 2 "" 1 Q(2)' ' (3.6) 

a t = at, A = &*, r < 1. 

Then, from (E3D and (12^01) . we see that, for x G C/Sf , 

u f / <5r log(x — x* 



\og(x-x*) = 
n \ log or 

p (s) log(x — s)ds = O 

x*—a t 



log St 

It then follows from corollary 12. II and the properties (12.411) that the function g t (x) satisfies 
the following 

Proposition 3.1. For sufficiently small 5t, there exist 5 > such that [x* — a t ,x* + a t ] C 
13$ and that 

gl(x) +gt(x)-^- l -= v t (x) (JjL^ , x G [a u f3 t ]/ {Bf U B*) 

&( x ) +g i(x)-^- l -<0, xeR/ Aj^uK/3 t ]j, (3 - 7) 

<7*(a;) = log a; + 0(1), x — > oo. 
where I is the following constant 

f /-(^, t > 1; 
l « t , t<l. 

T/ie function v t (x) remains uniformly bounded in [a t , /3 t ] as St — > 0. TTiai is, «/ 

lim f+(x) = -u(x), 

r/ien -u(x) is uniformly bounded in [at, fit] . In particular, Vt(x) is zero when t < 1. 

The function g t (x) is analytic on C/(— oo, x*) and has the following jump discontinuities 

on (— oo, a*) and ((3 t ,x*). 



' ' - 1 °- : xe(-oo,o i ) ( 3 - 8 ) 



# + (x) - g_(x) = 27rz, i6(-oo,a t ) 

- ^(x) = 27TZ — , X e (fy, X*] 



n 
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we can now transform the Riemann-Hilbert problem with the function g t (x). 
Let T(x) be the following function 



Tlx) = e^Y(x)e~ n9 e n» 



where 03 is the Pauli matrix 



o"3 



1 
-II ' 



Then T(x) is a solution to the following Riemann-Hilbert problem 

1. T(x) is analytic in C/R; 

2. T+(x) = T_(x)J T (x), x e R; 

3. T(x) = I + 0(x~ 1 ), x-^00. 

where Jt{x) is the following matrix on R. 

e -n(g\(x)-gL(x)) e n{g t + {x)+g^(x)-V t (x)-\) 



Jr x) 







^n(g t + (x)-g t _(x)) 



X G 



(3.9) 



3.2.2 Opening of the lens 

We now perform a standard technique in the steepest decent method [3], [15], [16]. First 
note that, from (13. 7p . we see that Jt{x) becomes the following on the interval [a t , fit]- 



Mx) 



e -n(g t + {x)-g t _{x)) e 2D n {x) 



,n($\{x)-gt{x)) I ' 



x e [a t ,f3 t 



where D n (x) is the function 



n5t 



D n (x) = v t (x)— -, t > 1, D n {x) = 0, t < 1 
log or 



(3.10) 



which is bounded on [0^,$] under the double scaling limit (ll.7al) . 

Then from (13.81) . the jump matrix Jt{x) has the following factorization on [a f ,/3 t ]. 



Mx) 



-'' >-"> "•' I I -W-W I I n(v t (x)-2g t + (x)+l)~2D n (x) 



n(Vt(x)-2gUx)+\)-2D n (x) ^ I I e -2D n (x) Q 



As in [3], [15], [16], we can open a lens around the interval [at, fit] as shown in Figured] 
and define the matrix S(x) to be the following 

r Tlx), 



x outside the lens; 



S(x) = < 



T(x) 



T(x) 



1 0' 



— e 



nf V t (x)-2 ff '( a; )+ij-2 J D n ( a: ) 



(3.11) 







e n(y t (x)-2g t ( a; )+{)-2 J D„(aO j j ' 



x e L- 
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Then the function S(x) will satisfy the following Riemann-Hilbert problem. 

1. S(x) is analytic in C/R; 

2. S+(x) = S-(x)J s (x), xeR; (3.12) 

3. S'(x) = / + 0(x _1 ), x^oo. 

where Js{x) is now defined by the following 

I 1 °\ 



e 2D ™ (a;) \ 
W=( _ e -2D n ( a: ) o )> xe{a t ,p t ) 



(3.13) 



A n(2 fl *(*)-Vt(a)-f)+2D„(*)\ 

J 5 (x) = I e ^ I, xe{-oo,a t ) 

/ e -27rm* e n(2 fl *(a;)- V t (x) -l t ) +2D„(x) \ 

Js(x)=\ Q gW I, x6(/3 t ,x*) 

Then from ( 13. 7ft . we see that for some large enough n and t close to 1 such that St^^ = o(l), 
we have 

1. D n (x) is uniformly bounded on [at,/3t]. 

2. e «(VW-^(x)-|) ^ Q Qn R/ 5 p u ^ A] j 

3. e «(^(-)-^W+l) o on (d U C 2 ) / (Bf U 5 2 ) . 
Therefore, in the double scaling limit, the jump matrix Js{x) behaves as 

j s (x) - /, (x e R u d u c 2 ) I f(J Bp u K, A] j , 

= _ e -!D.(«) J ' X G 
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3.3 Parametrix outside of the points at, f3 t and x* 

We will now construct the parametrix outside of the singular points. We would like to find 
a solution to the following Riemann-Hilbert problem. 

1. S°°(x) is analytic in C/ (R U Bf ); 

2. S™(x)=S™(x)J°°(x), xGR; (3.14) 

3. S°°(x) = I + 0(x- 1 ), x^oo. 

where J°°(x) is the following matrix- valued function. 

/ -27riw* q \ 

J°°(x)={ e wJ, *e (&,*•) 

J°°(x) — [_ e -2D n ( x ) o )> xe(a t ,Pt). 

We should construct several scalar functions and use them to 'dress' the 1-cut parametrix 
constructed in [T3] so that it satisfies the Riemann-Hilbert problem f)3.14p with the jumps 

3.3.1 The limiting Abelian differential 

The discussions in sections 13.3.11 and 13.3.21 are only relevant when t > 1. Let us first 
construct a function F(x) with the following jump discontinuities 

F+(x) = -F_(x), xe[a t ,p t ] 
F+(x) = + 2vri, x G [A, x*] 

We have the following 

Lemma 3.1. The function F(x) defined by 



F(x) J f *•-<*)(*•- l08 V^lY^ (3.16) 



^/(s - a t )(s - /3 t ){s - x*) 



St* 


-a t 


X* 


-fit 


X* 


-at 


X* 


-fit 



+ 



x-Pt 



satisfies the following scalar Riemann-Hilbert problem. 

1. F(x) is analytic in C/[att, x*]; 

2. F+{x) = -F-(x), xe[a t ,/3t\; 

3. F + (x) = F_(x) + 2?rz, x e [A, x*}; 

4. F(x) = \og(x - x*) + 0(1), x^x*; 

5. F(x) = Fqj+O^ 1 ), x^oo; (3.17) 
- 1 

^0 = Wg ; 

1 



X* 


-a t 


X* 


-fit 


X* 


-a t 


t* 


-fit 



6. F(x) is bounded as x approaches a t or (3 t 
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where the square roots and logarithm are chosen such that the branch cut is on the negative 
real axis. Also, to avoid ambiguity, we set 



x 



a t _ y/{x- a t )(x - p t ) 



x-& (x-Pt) 2 

Proof. First let us show that the only singularity of F(x) is at x*. The function F(x) 
can become singular when the argument becomes zero or infinity. The numerator and the 
denominator of the argument can become zero if 



x* — a t \x — a* 

± \ i- = 



x* — Pt \ x — p i 
which implies 

x* — a t x — a t 

x* - p t x - Pt 

X = X* . 

Since the denominator is non-zero at x = x*, we see that the denominator does not vanish 
for all x G C. Near x = x*, we can expand the numerator and the denominator in a power 
series of x — x*. 



X* 


-a t 


.('* 


-A 


' X " 


-at 


X* 


-Pt 



fx- 


-at 


X 


-fit 


'x- 





c (x - x*) + 0((x - x*y), 



W x- I 

where Cq is the following constant 



d I x—at 

dx\\J x-fkj \ x=x * = p t - at , 

C °~ ~ 4(x* - p t )(x* - a t ) ^ 

x*—fit 

Therefore near x = x*, F(x) behaves like 

F(x) = log (a; — x*) + 0(1), x — > x*, 

this proves 4. 

The other points where F(x) can be singular are the points x = at, Pt or x = oo, where 



the argument may become infinite. However, near x = at, the function y^z^ remains 

finite and therefore F(x) does not have singularity near it. Let the argument inside the 
logarithm of ( 13. 161) be 



X * 


-a t 


X* 


-fit 


X* 


-a t 


x" 


-fit 



X- 


-a t 


x- 


-fit 


x- 


-at 


x- 


-fit 



= \ ^ "-^=. (3.18) 

+ 
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Then near x = ft, $(x) behaves like 



$(x) 



which is bounded and non-zero as x — > /3 t . 

We now consider the point x = oo. Near x = oo, we can rewrite <&(#) as 



' (a:*— Qj) 
x*-/3t 



0*-nt) , Z 1 -^ 
x*-/3t ^ y 

Therefore the asymptotic behavior of 3>(x) near x = oo is given by 



ti- 


21 






ft 


X 


ll- 









-at) 


X* 


~Pt 




-a t ) 



,x*-q t ) I 1 

x*-/? t 1 

Hence has a singularity at x = x* only and this proves 5 and 6. 

We will now study the jump discontinuities of F(x). First note that F(x) can only 
have jump discontinuities outside [o; t ,/3 t ] if 



X* 


-at 


X* 


-fit 


>x* 


-a t 


X* 


~Pt 



x-at 
x-Pt 



$(ar) = V e R. 

, / X-Q« 

Simple calculations shows that x G K. Therefore -F(rr) can only have jumps on the real 
axis. 



We will first consider the jump on [a t ,(3 t ]. On [ct t ,/3 t ] the square function y fz^ has 
the following jump discontinuity. 



x — cxt\ I x — at 



x - ft J \ V x - ft 



Hence F-(x) is given by 



'x* 


-a t 


X* 


-ft 


'x* 


-a t 


X* 


-Pt 



X*-Qf 

x*-ft 



= log = log = -F+(ar). 

, / / x-q t \ / x*-qj _ 

This proves property 2. 



x- 


-a t 


x- 


-Pt 


x- 


-a t 


x- 


-Pt 
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We now consider the jump on [/? t , sc*]. The function in (13.181) is real on M./[ott, fit] 
and we need to show that it is negative on (fi t ,x*) and positive elsewhere. Let us first 
consider a point x € (fit, °o). The the denominator of is positive for x G oo). To 
study the signs of the numerator, let us consider its derivative 



d ( j x* - a t fx - a t \ 1 / x - fi t fi t - a t , . 

we see that the numerator is a strictly increasing function on (fit, oo). Since it vanishes at 
x = x*, we see that it is negative on (fi t ,x*) and positive on (x*, oo). 

Now let us look at the sign of $(x) on (— oo, at). First note that, on (fi t , oo), the square 



root J ^rq£ is positive and strictly decreasing and hence it is greater than 1 on (fi t} oo). In 
particular, we have 



x* — fi } 



Now let x £ (— oo,Q! t ). In this region, the square root yf^f^ ^ s positive and strictly 

decreasing from (13.191) . Near — oo, it approaches 1 while at a t , it becomes zero. Therefore 
on (— <x>,at), it takes values between and 1. 
Therefore we have 





- at 


X* 


-A 


'x* 


- at 


X* 


-fit 



' X 


- a t 


X 


-A 


1 X 


- a t 



> 1 — 1 = 0, x E (— oo, a t ) 



+ \ w > 0, xe(-oo,a t ). 

V x — pt 

Hence the $(x) is positive on (—oo,a t ). Summarizing, we have 

$(x) > 0, x E (-oo, a t ) U (x*, oo) 
$(x) < 0, x G {(3 t ,x*). 

This proves 3. 

Since F(x) cannot have any jump discontinuities and singularities other than the ones 
that are considered here, property 1. is true. □ 

Remark 3.1. The function F(x) can be thought of as the limit of an Abelian integral on 
an elliptic curve. Let us consider the following elliptic curve 

z 2 = (x- a t ) (a: - (3 t ){(x - x*) 2 - a 2 ) (3.20) 

and define the a and b-cycles of this curve as in Figure® Then the normalized holomorphic 
Abelian differential on this curve is given by 

Q(x) = . ° dX = (3.21) 

y/(x- a t ){x - f3 t )((x - x*) 2 - a 2 ) 
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Figure 2: The a and b cycle of the elliptic curve (13.201) . 



for some constant C such that 



Q{x) = 1 



In the limit at — > 0, the curve becomes degenerate and the Abelian integral J x fl(s) degen- 
erates into the function F(x). 



3.3.2 Scalar function with jump D n (x) 

We will now seek a scalar function K n (x) that is bounded at infinity and has jump 2D n (x) 
on [a t ,/3J. (cf. the Szego function used in 1'23]. [32] ) 

We shall construct a function K(x) that satisfies the following Riemann-Hilbert prob- 
lem. 

1. K(x) is analytic on C/[a t , /3 t ]; 

2. K+(x) = -K^{x) + 2D n (x), xe[a t ,Pt]l 



3. K{x) =K + 0(x 



-In 



x — > oo, 



Kn 



1 



2m 



at 



2D n {s)ds 



yj{s- a t ){s - A) 



This function can be constructed by the use of Cauchy transform easily 

2D n (s)ds 



K{x) 



^(x-a t )(x-[3 t ) f 



2m 



at V ( s _ a t)( s ~ A) (s-x) 



(3.22) 



(3.23) 



Lemma 3.2. The function K{x) defined in \3.23\l satisfies the Riemann-Hilbert problem 

m 



Proof. Let C(f) be the Cauchy transform 



C(/)(x) 



f(s)ds 
2m J at s-x 
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Then from the Plemelj formula (See, e.g. [26]), we have 

i r Pt f{s)d s 



C=(f)U) :- 2 .f ~ I .rE\n,. i n . 



2ni ./ s - x 



where the principal value is taken in the integral on the right hand side. Taking into 
account the change of sign of y/(x — a t )(x — (3 t ) across [a t , (3 t ], we have 

K + (x) + K_(x) = (y/(x - a t )(x - ft)) ? ; 2Dn(x) 
= 2D n {x). 

To see that K{x) has the desired property at x = oo, let us write the factor in a power 
series near x = oo. 

— = ~(l + - + 0(x- 2 )). 

s — X X \ X / 

Therefore the function K(x) behaves as follows as x — ■> oo. 

v/(x-o 4 )(x-/3 t ) ( 2D n (s)ds ^ 



2nix 



V 



'<* (V(s-a t )(s-/3 t )) 



0(x 



2D n {s)ds 



2mJ at (y (s _ a<)(s _ A) 
which is property 3. in (13.221) . 



+ o(x- 1 ; 



□ 



3.3.3 Parametrix outside of special points 

We are now in a position to construct the parametrix outside of the special points. First 
let us consider the following matrix. 

f(x)+j(x)- 1 ■y(x)--y(x)- 1 

Tl(x) = I ..^ 2 .,^-! ..^, 2 .^- 1 1 (3.24) 



2 _ 2i 
f(x)— 7(a) 1 7(x)+7(a) 1 

2i 2 



where 7(2) = jt(%) is defined by 



7(x) 



x- p t 



X — OL t , 

Recall that this matrix satisfies the following Riemann-Hilbert problem [T5], [T2], [3]. 

1. n(x) is analytic on C/[a t , /3 t ]; 

i N 



2. n + (x) = n_ x 



-1 o 



, x G [a t , A 



3. n(x) = l + 0{x~ x ), x -> oo. 



Birth of new cut double scaling limit 



29 



We can now combine U(x), K(x) and F(x) to form our parametrix outside the special 
points. The main result is the following. 

Proposition 3.2. The matrix S°°(x) defined by 

S°°(x) = e ( x o + ( ut ~"') F o) CT 3TT^ x ) e -(-^( :1 ')+(« t -« t )-p , (^))o-3 (3 25) 



satisfies the Riemann-Hilbert problem jjj3.14\ ) an d A3.15\) . In particular, when t < 1, both 



u l — u l and K(x) vanishes and the above equation reduces to 

S°°(x) = U(x), t < 1. (3.26) 
Proof. First let us consider the asymptotic behavior near x = oo. We have 

S°°( X ) = e (Ko+(u*-u*)F )* 3 + Qix- 1 )) e ""(^ + («)^3 

= (/ + 0(x- 1 )) 

This proves 3. in (I3.14j) . Next we consider the jump discontinuities. It is given by 



S T{x) = S™(x)( eAl ^ JL), xG(-oo,x*)/K,A] 



e" 

i e A 2(x)\ 
S?(x) = S™(x)l^_ e y Mx) e Q J, xe(a t ,(3 h 

where Ai(x) is given by 

A 1 (x) = (u* - u 1 ) (-F + (x) + , x e x*), 

A 2 (x) = (tZ* - u*)(F+(x) + F_(z)) + K+(x) + K_{x), xe (a t ,p t ). 

From (13.171) and (I3.22p . we see that Ai(x) and A 2 (x) are in fact the following 

At(x) = -lixiiu 1 - if), x G (Pt,x*) 
A 2 (x) = 2D n (x), xe(at,Pt). 

Since u* is an integer, we see that 



(3.27) 



x G {(3 t ,x*). 



Substituting these back into ( 13. 27ft . we see that the matrix S°°(x) does indeed satisfy the 
jump conditions (13.151) . □ 

Remark 3.2. The appearance of the degenerate Abelian integral F(x) $3.16}) in the parametrix 
$3. 25\) comes from the appearance of the elliptic theta function in the 1-cut parametrix. Let 
the 1-cut parametrix outside of the special points be M°°(x). Then M°°(x) is given by [It 







f 7+7" 


- 1 e(w(x)-(u t -u t )+d) 


7-7-1 e(-W(x)-(u t -u t )+d)\ 


M°°(x) 




. 7-7" 


6{W(x)+d) 

- 1 e{w{x)-{u t -u t )-d) 


-2i 0(-W(x)+d) | 
7+7- 1 e{-W{x)-{u t -u t )-d) 1 






\ 2i 


e(w(x)-d) 


2 d(-W(x)+d) J 






( 7-7" 


- 1 9(-W(x)-(u t -u t )+d) 


7+7- 1 e(w(x)-{u t -u t )+d) 


M°°(x) 




2i 
7+7" 


6{-W{x)+d) 
- 1 B(-W(x)-(u i -l^)-d) 


2 8(W(x)+d) 

7-7- 1 e(w(x)-(u t -u t )-d) 



3(z) > 



Q{z) < 0. 



9{W(x)+d) 2i 8(W(x)-d) 
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for some scalar constant d and constant diagonal matrix H , where W(x) is the Abelian in- 
tegral J x Q(s) $3.21\) and 9(z) is the elliptic theta function. One can think of the parametrix 
S°°(x) as a degenerate version of M°°(x) as the branch cut [x* — a t ,x* + a t ] in $3.20\) is 
closing up and the curve degenerates into a genus zero curve. In this case, the Abelian 
integral degenerates into the function F(x) and while the theta function degenerates into 
an exponential function. In the multi-cut case, one could apply the analysis similar to those 
in 1^, J2iy to obtain degenerate hyper- elliptic theta functions and use them to construct 
the suitable parametrix. 

3.4 Parametrix near the edge points a t and (3 t 

At the edge points a t and (5 t the approximated density vanishes like a square root 

and the local parametrices S ±2 (x) near these points can be constructed by the use of Airy- 
functions. Such construction has been done many times in the literature and we should 
not repeat the details here. An interested reader can consult [T5], [TB] , [3] for example. 

3.5 Local parametrix near the critical point x* for t > 1 

We will now consider the parametrix near the critical point x* . As in [TH] , the parametrix 
will be constructed out of the monic orthogonal polynomial 7r^ t (() of degree vf and weight 
e - ^, where 2z/ is the order of vanishing of 2h{x) — V(x) + I at x* . 
We would like to construct a parametrix S x *(x) in Bf such that 

1 . S x * (x) is analytic in Bf / (Bf n R) ; 

2. S x + \x) = S x _\x)J s {x), x e Bf n R; (3.28) 

3. S x *(x) = (I + o(l)) S°°(x) as n -> oo, t -> 1, uniformly in dBf . 

3.5.1 Conformal map in Bf 

Let us define a conformal map ( = f(x) that maps the neighborhood B$ into the complex 
plane, such that, as n — >• oo, the boundary of Bf is mapped into infinity. We will define 
( as follows. 

C = f(x) = (-n (2h(x) - V(x) -l))^, (3.29) 

where the ^th-root is chosen such that the intervals [x* — 5, x*] and [x*, x* + 6} are mapped 
onto the negative and positive real axis respectively. This is possible because h(x) — — | 
vanishes to order 2v at x* and that it is real and negative on the interval [x* — 5, x* + 5] 
due to (g~n$ . 

Since h(x) — ^y^- — | vanishes to order 2v at x*, the function ( is of the form 

t = n&(x-x*)<p(x) (3.30) 
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such that <^(x) is independent on n and <p(x*) 7^ 0. By choosing 5t and S smaller if 
necessary, we can assume that (p(x) and hence ( is conformal inside the neighborhood B'f . 

Then ( maps the neighborhood Bf into the complex £-plane such that the boundary 
of Bf is mapped into infinity. 

Let us now define the constant Z t and function r t (x) by 



z, = n \ u\\r s — — — Yt) u \^ 0gtp { x *) + ^; lo s n 

n (2^(1) - - I) - 2Z t + C v ~ 2«« log (n&<p(x) , , , 

nix) = -± ^ S V- 61 -< 

Note that Tt(:r) does not have a pole at x = x* and that by taking St and S smaller if 
necessary, we can assume that n(x) is analytic inside Bf . 
Then it is easy to see that (, n(x) and Z t satisfy 

nfa-m-£)-JZ + &< + ^< + ± «3, 2 ) 
Moreover, we have the following 

Proposition 3.3. As n — > 00 under the scaling \l.7dty , the constant Z t and function n{x) 
are of order 

Z t = -—log logn + 0(1), 

7, N (3.33) 

n(x) = 

uniformly in Bf . 

Proof. The first part of the proposition follows from (12.431) . By (I2.43[) and ( 13 .3 1 f) . we have 

St 



Z t = nSt(f)(x*) — u l ^log<^(:r*) + — logn^ + O 



log St 



Then by using (13.41) to eliminate nSt, we obtain the following 



^ = (logn + logtft) - u'log^*) + O ) • (3-3 I ) 

Now by taking the logarithm of (ll.7al) . we see that, 

logn + log 5t = log logn + 0(1). 
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Hence we obtain 

u* 

Z t = -— loglogn + 0(1). 

This proves the 1st equation in (13.331) . 

Given the first equation in (13.331) , the second equation in (13.331) now follows easily from 
the definition of C (1^291) . r t (x) (EPTD and the relation fl2T42l) from corollary E2 □ 



3.5.2 Construction of the parametrix 

We should now construct the parametrix by using orthogonal polynomials with weight 
exp(-C 2 " + T t (x)C). 

Let 7r£(C r) be the monic orthogonal polynomial of degree k with respect to the weight 
' <(CK(C) exp {-C 2u + <) dC, = K(r)5 k3 , 



where h v k {r) is the normalization constant as a function of r. Let us denote by ty u (£,s) 
the following matrix constructed from the orthogonal polynomial ir^ t ((,T). 

( n u (C t) W < f (^)cxp(^ 2 -+r C ) \ 

<C ' " L l(TW l(f ,r) %gj,^"-^J ' <X35) 

\ W — 1 ^ / IT — 1"' / 27TI JM S — X / 

where = -pf^y. 

Then the matrix f((,r) satisfies the following Riemann-Hilbert problem. 

1. *"(C, t) is analytic on C/R; 

2. n( C,r) = *-(C,r)(j «P("<* + CeRi 



3. M»"(C,t)=(/ + 0(C- 1 ))( C °' ( °,.), C 



oo. 



Let .E'(x) be the following matrix- valued function, 

E(x) = S^ixX^'^e^ 3 (3.37) 

Then from (I3.25f) we see that 

E{x) = e ( Ko+iut -^ )F °) a3 U{x)e- K ^ a3 ^-(u*-*)^)^^-*)^ ( 3 . 38 ) 
From property 4. of (13.171) . we see that the factor 

e -(«*— u*)F(a;)cr3^-(u*— u')(T3 

is analytic inside Bf. Then, since both K(x) and H(x) are analytic inside Bf* , the function 
E(x) is analytic inside B% . Hence we have 
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Proposition 3.4. Let the matrix S x (x) be 

S x *{x) = E{x)V u {(,T t (x))(- ut ' 73 e- Zt ' 73 , xeBf, (3.39) 
Then, under the double scaling limit fll. 7a\j , S x * (x) satisfies the conditions 



1. S x * (x) is analytic in Bf/ (Bf n R) ; 

2. Sf (x) = Sf (x)J s (x), x G Bf n R; 

3. .S" (.,-) = ( I + O I J, ) ) S' X V) «* /' - x. aniformhi in OBj 



(3.40) 



Proof. The properties 1. and 2. follows immediately from (I3.32p and property 2. of (13.361) . 
We should now prove property 3. 

At the boundary of Bf , we have ( — > oo and hence the function S x *(x) behaves as 



S x *(x) = S 0O (x)C (uW)CT3 e ZtCT3 (/ + 0(C -1 )) e - Zta % {ut - ut)a3 , ( 



oo. 



From (13.291) . we see that ( 1 = 0{n ^) at the boundary of Bf , hence the above equation 
becomes 



S x (x) = S^ix) [I + O 



(log 71 ) 2* 

l-2|«*-tt*| 

n 2" 



where the second equality follows from the fact that S°°(x) is bounded in as n — ► oo. 
This proves the proposition. □ 

3.6 Local parametrix near the critical point x* for t < 1 

We will now construct the parametrix in Bf when £ < 1. In this case, the parametrix can 
be constructed from the Cauchy transform |16j . 

3.6.1 Conformal map in Bf 

We will use the same conformal map (13.291) defined in section I3.5.1L However, the function 
r t (x) and the constant Z t are now defined to be 

_ Y^l _ k 

2 2 

7 (3 4r 

. n (2g t (x) - V t (x) - l t ) - 2Z t + ( 2u 
n(x) - - . 



Z t = n( g (x 
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Then by ( 11.7bl) and the Buyarov-Rakhmanov equation (11. lip , we see that Z t is of order 

Z t = 0(n8t) = 0(n 1 - k ) (3.42) 



As we see in ( 11. 17b II . the limiting kernel will be of order e Zt . However, one should bear in 
mind that Z t is negative and therefore the term e Zt is bounded as n — > oo. 
We can now deduce the order of T t (x) from (II. lip . (11.7bl) and (I3.42p . 

(3.43) 

[X — x*)tp(x) 

From the definition (13.41 j) . we see that (, Z t and Tt(x) together satisfies the following 

V t (x) l,\ _ C 2 " , T,(X)( 



3.6.2 Construction of the parametrix 

Let us now construct the local parametrix by using Cauchy transform (cf. PS]). Let 
ty((,T t (x)) be the following matrix. 

A 1 f cxp(-C 2 "+r t (x)C) \ 

ty((,T t (x)) = f J 2^ Jr ^ ds \. (3.45) 

Then \I/(£, r t (x)) is the unique solution to the following Riemann-Hilbert problem. 

1. r t (x)) is analytic on C/R; 

2. <MC,n(*)) = MC,n(*))(j ex P(-^ + <)^ CgM . (346) 

3. y(C,T t (x)) = (l + 0(C 1 )), C^oo. 

We can use ty((,Tt(x)) to construct the local parametrix S x (x). Let us define the matrix 
E(x) to be 

E(x) = S co {x)e Zta3 (3.47) 

then E(x) is analytic and invertible inside of Bf . From (I3.44p and (13.461) . we have the 
following. 

Proposition 3.5. Let the matrix S x (x) be 

S x *(x) = E{x)^((,T t (x))e- ZtCT3 , x G Bf, (3.48) 



Then, under the double scaling limit (1.7b), S x *(x) satisfies the conditions 

1. S x * (x) is analytic in Bf/ (Bf n R) ; 

2. Sf(x) = Sf(x)J s (x), x G Bf fl R; (3.49) 

3. ^(x) = + O (n~^ j \ S°°(x) as n — > oo ; uniformly in dBf . 
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Figure 3: The contour £ on which R(x) is not analytic. 

Proof. As in the proof of proposition 13.41 properties 1. and 2. are clear from (I3.44p and 
(13.461) . Let us take a look at the condition at the boundary of B$ . We need to be careful 
as Z t may contain powers of n in it. At the boundary of Bf* , we have ( — > oo, and 



S x * ( x ) = 



1 0{e 2Zt C l ) 
1 



but from the expression of Z t (13.411) and the property of the equilibrium measure (12. 3p , we 
see that Z t < and hence e 2Zt is bounded. Hence we have 



yoo (x) 



S x \ x ) = S°°{x) (7 + (™^)) = (j + S c 
This completes the proof of the proposition. □ 

3.7 Last transformation of the Riemann-Hilbert problem 

Let us now define R(x) to be the following function. 



S(x) (S n (x)) 1 , x inside B 



n. 

R ^ = { S(x) {S°°\x)y 1 \ x outside of B r 5 \ (3 ' 50) 

where r\ = —2, r2 = 2, r% = x* and S ±2 (x) are the local parametrices near the edge points 
at and fit- Then the function R(x) has jump discontinuities on the contour £ shown in 
Figure [3J 

In particular, R(x) satisfies the Riemann-Hilbert problem 

1. R(x) is analytic on C/S 

2. R + (x) = R_(x) J R (x) 

3. R(x) = I + 0(x- 1 ) 1 x-*og [ ' ' 

4. R(x) is bounded. 

From the definition of R(x) ( 13.501) . it is easy to see that the jumps Jr(x) has the following 
order of magnitude. 

' I + Oin- 1 ), x G dBg 2 UB 2 s ; 

/ 2\u t -u t \-l u '\ * 

l + 0[n (logn)^ , x G dBf , t > 1; 

J*(*) = > x\ ' ' ( 3 - 52 ) 

J + 0(n-^), xedB%, t<l; 



I + O (e ™ 7 ) , for some fixed 7 > on the rest of S. 
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Since \u l — u t \ < |, for sufficiently large n, n 1 2^ 1 (logn)^ and n ^ are small. Then 
by the standard theory, p2], [15], [16], we have 

i?(x) = / + [n ~v (logn)^ , t>l 

V / (3.53) 

= / + O (n~^ V t<l. 

uniformly in C. 

In particular, the solution S(x) of the Riemann-Hilbert problem (13 .12ft can be approx- 
imated by S°°(x) and S n (x) as 

cm - S R(x)S ri (x), x inside B r /; 

* [X> ~ \ R{x)S°°{x), x outside of B r s \ {6 ^ > 

When t > 1, this approximation becomes poor as \u l — w*| gets close to |. However, if we 
restrict our attention to a small neighborhood of x* such that 

z = (x — x*) n^(p(x*), (3.55) 

is finite, then we can still use this approximation to obtain the asymptotic kernel (11.1 7aj) . 



4 Asymptotics of the correlation kernel 

We should now compute the kernel using the the asymptotics obtained in section [3j Recall 
that the kernel and the solution Y(x) of the Riemann-Hilbert problem (13. lj) are related by 

4.1 Asymptotics of the kernel when t > 1 

First let us recover the asymptotics of Y(x) from that of S(x). By reversing the series of 
transformations (13.111) and (13.91) . we find that, for x G Bf , the matrix S(x) and Y(x) are 
related by 

Y{x) = e n ^ 3 S(x)e n (. 9t{x) -^h, x G Bf. 
We now use the estimate (13.541) and the expression of (13.39!) to obtain 

Y(x) = e^ CT3 J R(x) J E;(x)vI/-(x,n(x))e n ( 9t(x) -*- utlogC - Zt ) CT;i , x e B f. 
Now from (13.321) . we see that the above is equal to 

Y(x) = e n ^ 3 R(x)E(x)^ u (x 1 r t (x))e( n ^"^ +I ^ K )' J3 , x G Bf. (4.1) 
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Let us now study the behavior of E(x) and R(x) in the vicinity of x* when z defined by 
(13.551) is finite. First let us consider E(x). From ( 13.381) . we see that E(x) is analytic inside 
the neighborhood Bf . Then from the power series expansion of E(x) inside B$ , we obtain 

E(x) = (e° + E 1 zn~^ + O (n'^yj n^ a3 (logn)~^ CT3 , (4.2) 

for some constants E° and E 1 that are bounded as n — ► oo. 
Now consider R(x). Let m be the biggest integer such that 



m 



(l - 2\u l < 1. 



Then from the Riemann-Hilbert problem (13. 50f) . we see that Rix) is analytic inside -Bf*, 
hence in terms of z, we have the following estimate 

/ 2|n t -iT t |-l 2it\' > 

R(x) = I + > [ n 2 " (log") 2 " J 

i=i V ' (4.3) 

+ O {^zn^^ (log n)^ J + O fn - ^ J , 

( 1 U* \ 1 ll 4 

where O ( zn~^\\.ogn)^ J denotes z dependent terms that are of order n~2U (logn)^ . The 

constants Aj are finite in the limit n — > oo. 

In particular, from (14.21) we see that E~ 1 {x')E{x) satisfies the following 



z' — z 



E (x')E(x) = n — — a3 (logn)^ 3 \ I + O ^ — —J J n^^ CT3 (logn)~^ ff3 , (4.4) 

while R~ 1 (x')R(x) satisfies the following estimate 

R-^x^Rix) = 1 + [{z' - z)rr^(logn)^) (4.5) 
Hence the product E~ 1 (x')R~ 1 (x')R(x)E(x) satisfies the following estimate 

E- 1 {x')R- 1 {x')R(x)E{x) = 1 + \{z - z)n ^ (logra) ^ ) (4.6) 



If we now substitute (14.11) and (14.61) back to (13.21) . we obtain the following estimate for the 
kernel 
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where V(() = -( 2u + T t (x)(. 

Now recall that from (13.331) and (13.301) . we have 

lim ( = z, lim r t {x) = 

n-^oo n— >oo 

If we now take these into account and substitute (13.35}) into (14.71) . then we obtain the limit 
of the kernel as 

hm —L^K n x<) = < gWfe(f) - m-^) 

n,N^^ x *) n l n ' NK ! u - 1 2m{z-z>) (4.8) 

where u = lim^oo u* and = «^(0). 

To complete the proof of the first part of theorem 11.11 we need to show that (p(x*) = 



Q(x')y/(x*)*-4, 
2v 



i 

21/ 



. This can be seen from the expression (12.51) . From (12. 5p . we have 



/—, — r f p(y) , V'(x) , V'(x) 

Vo(x)= / ^ L dy+—^ = -h'{x) + — \-L, 

then from the fact that 2h(x*) — V(x*) — I = and the expressions of of q(x) ( \2.6h and ( 
(13.291) . we see that, upon integration, we have 

/ *\2u 2v I *\ / *\2i/Q( X )v( a '*) 2 — ^ 

n(x — x ) ip (x ) = n(x — x ) 

this completes the proof of theorem II. II for the case t > 1. 

4.2 Asymptotics of the kernel when t < 1 

We will now use the local parametrix S x (x) (13.481) constructed for t < 1 to compute the 
kernel. In this case, the solution Y(x) to (13.11) is given by (14. ip with ^/ u (x,r t (x)) replaced 
by ^/(x, T t {x)) and r t (x) defined by (I3.4ip . 

Let z be the variable defined by (I3.55P and assume that z is finite. Then by using the 
power series expansion of E(x) and R(x) inside Bf, we obtain 

E(x) = I U(x*) + U'(x*) Z —^ + O (zn-^) ) e Zta \ 

\ ip(x*)n^ V / J (4.9) 

R(x) = I + R°n-^ + O (n~^\ + O (zn~^\ 

where we have used (13.261) to replace S°°(x) by U(x) and O (^zn~^^ denotes z dependent 
terms with order n~^. Therefore the product E~ 1 (x')R(x')R(x)E(x) is of order 

E-\x')R{x')R{x)E{x) = e~ z ^ ( I + n ~V)n(z*)0z Z <> + Q (l^j£\ ) 

y ip(x*)ri2v V nzv J J 

(4.10) 
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From (13.451) . one can easily check that 



(4.11) 



o J V o ) ' 

where V(£) = ~C 2u + r t (x)C- 

We then substitute (14. 10p and (13.451) into (13. 2p and arrive at 



p ( 2Z * + — + —2-J /l\ , , . xx 

^(x,x') = — (o i)n- x (x*)n(x*) (J) (1 + 0(71--)) 



2m 

This gives the double scaling limit of the kernel 

hm e iZt K rhN (x,x') = e ? -- — - — . 

n,N-* 00 87T \X* — p t X* — a t J 

This completes the proof of theorem 11.11 
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